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Abstract. - 

Motivated by recent experiments, in which the Kondo effect has been observed for the 
first time in a double quantum-dot structure, we study electron transport through a system 
consisting of two ultrasmall, capacitively-coupled dots with large level spacing and charging 
energy. Due to strong interdot Coulomb correlations, the Kondo effect has two possible sources, 
the spin and orbital degeneracies, and it is maximized when both occur simultaneously. The 
large number of tunable parameters allows a range of manipulations of the Kondo physics - in 
particular, the Kondo effect in each dot is sensitive to changes in the state of the other dot. For 
a thorough account of the system dynamics, the linear and nonlinear conductance is calculated 
in perturbative and non-perturbative approaches. In addition, the temperature dependence of 
the resonant peak heights is evaluated in the framework of a renormalization group analysis. 



Introduction. Recently there has been substantial interest in many-body and correlation 
effects in ultrasmall semiconductor quantum dots. The dots may have strong electron-electron 
interactions, characterized by the charging energy Eq, and large level spacing S (typically 
Ec > 6) At low temperature and strong coupling T between the dot and the leads, 
fcer <C r, quantum fluctuations of the charge and spin degrees of freedom strongly affect the 
transport through the dot |^|-^. The spin fluctuations lead to the Kondo effect, which has 
been verified in experiments on single quantum dots 0. 

The Kondo effect has recently been discovered also in a double quantum-dot structure |^ . 
Motivated by these experiments, we consider in the present work a system of two capacitively- 
coupled quantum dots, depicted in Fig. |l|a). The level spacings in each dot are large, and 
effectively just one level per dot is coupled to the two reservoirs. For strong tunnel coupling 
and spin-degenerate levels, each dot would separately display the usual spin-Kondo effect in its 
I — V characteristics. In the present case, the interdot interaction is assumed strong rendering 
the charge states ni and n2 of the two dots strongly correlated. Consequently, the Kondo 
effect of each dot is sensitive to and can be manipulated by voltages applied to the other dot. 
The strong interdot interaction has also another, more dramatic effect: the Kondo effect can 
arise from the orbital degeneracy (the two dot levels tuned to resonance) even in absence of 
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Fig. 1 - a) Schematic of the capacitively-coupled double quantum dot. The bowed arrows denote 
possibility for tunnelling; there is no direct coupling between the dots, b) Part of the plane spanned 
by the gate voltages - the charge configurations ni , 712 minimize the charging energy within the 
hexagonal areas. In the Coulomb blockade regime, electron transport is enabled through dot 1 (2) 
along the thick (dashed) lines, while, in the rest of the plane including the thin solid lines, charge 
transfer is only possible via cotunnelling and other higher-order processes. 

the real spin degree of freedom The present model, with two spin-degenerate levels, allows 
us to study the interplay of the two Kondo effects as a function of different voltages as well as 
the level splitting. In absence of orbital degeneracy, the differential conductance can exhibit 
satellite resonances for finite bias voltages Vi ; these may be easier to observe than split peaks 
in the case of single dots {e.g. due to Zeeman splitting j3[Q or a second level in the dot ||^). 

The paper is organized as follows. First the model for the double-dot system is introduced. 
Then we express the current and conductance in terms of the spectral density of the dot 
states in the framework of a real-time diagrammatic technique. We calculate the linear and 
nonlinear conductance through each dot as a function of the two transport and the two gate 
voltages. This calculation is complemented by a poor man's scaling analysis of the temperature 
dependence of the linear conductance. 

The model. The double-dot system shown in Fig. ^a) is governed by the Hamiltonian 
H = Hrcs+Hdot+WT, where Tires = J2kricr '^kricr^kria'^kricr dcscribcs uou-interactiug electrons 
in the reservoir ri {r — L,K and i — 1,2), Tidot — J2icr ^ia-'^la'^ia- + Uni,n2 refers to the dot 
electrons, and = J^kriai'^rf '^kricr'^icr + H.c.) is the tunnelling Hamiltonian. For each 
barrier we introduce the tunneUing strength Tri — 27r^^ \TI^^\'^S{uj — Ekria) ~ 27r|Trippo, 
which is assumed to be independent of energy in the range of interest (po is the density of states 
in the leads). The electron-electron interaction, for electrons in the dot i, is accounted for 
by the charging energy Uni,n2- In general, the interaction depends on all the voltages applied 
to the system, but for symmetrically applied transport voltages it is determined by the gate 
voltages Vgi and Vg2 alone. 

In the plane spanned by these gate voltages, within the hexagonal areas depicted in Fig. |^6) , 
the ground state is reached for the charge configurations (71,1,712) indicated in the figure, and 
electron tunnelling is suppressed by Coulomb blockade effects. In the figure, there are three 
kinds of boundaries in the honeycomb lattice. Electron transport through dot 1 (2) is only 
possible along the boundaries indicated with the thick (dashed) lines, while the thin lines 
correspond to degeneracy between the two dots, in particular, ^lo- = ^2(t along the a-b line. 
As in the experiment, where the two dots lie close on top of each other, both the intradot and 
interdot interactions are assumed to be strong. Hence, for low temperature and bias voltages, 
k^T, eVi <^ Ec, it suffices to consider three charge states. Below we focus on the circled region 
in Fig. |l|6), i.e. the states (0,0), (0,1), and (1,0). In this restricted basis, the states can be 
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newly defined to include the charging energy C/„i,„2, e.g., for an electron with spin a in dot 
1, Cicr + Uifi — Uofi £icr- We observe that the resulting Hamiltonian H has a one-to-one 
correspondence to a single-dot model with two (spin-degenerate) orbital states, distinguished 
by an orbital quantum number i. The index i corresponds to the spatially separated dots and 
is therefore conserved in the tunnelling processes. 

Transport. We aim to evaluate the differential conductance Gi{Vi,V2) = dh/dVi [G = 
Gi +G2) from the current through the dot i. The dc-current can be expressed as [|| 

T T 1'^ 

= f-^^^ / Mfmi^) - fUioU^Acj). (1) 



Li -I- 1 Ri 



Here frii'^j) are the Fermi functions of the leads with chemical potentials /^ri- The spectral 
density Ai^{uj) = [C^{ui) — G^{uj)]/2TTi of the state ia corresponds to a local density of states 
in the quantum dot(s). Its shape is reflected in the differential conductance. 

In the weak-coupling limit, the lowest-order perturbation theory yields (to) oc d{uj—eic^) 
which describes the suppression of the current, i.e., Coulomb blockade effects. In second order, 
or 'cotunnelling' regime, the spectral density takes the form 

^^^^^('"^ = ^""JI^^^T-Tm^ {^^°"' +p!°^)(7+(^) + 7r(^)) (2) 

where are the occupation probabilities for the states ia as obtained in the lowest-order 
perturbation theory ||l2| and 271^^(10) = J^r^rif^ii^^)- As the prefactor in Eq. (||) suggests, 
j4jcr(w) only decays algebraically as l/(a; — eio-)^ leading to a non-vanishing conductance 
everywhere in the plane in Fig. The resulting linear conductance - denoted as g in the 
following - as a function of Ae = £20- —^lo- shows a maximum for Ae = 0, i.e., on the line a-b in 
Fig. |l|&) (this straight-forward result is not displayed). As a function of the transport voltage 
V, for symmetric applied voltages V = Vi = Vjn — Vi,i for both i, the differential conductance 
exhibits smooth variations, which are displayed as the two lowest curves in Fig. ^p). 

For stronger coupling, further higher-order processes become important. We account for 
them within the so-called resonant-tunnelling approximation Q| , which amounts to a resum- 
mation of a certain class of diagrams (with maximally two-fold off-diagonal density matrix) in 
all orders in the tunnelling amplitudes. To gain some insight into the results, let us consider 
again the symmetric applied voltages, equal barriers, Tri — F^, and degenerate levels Sia = £q 
{i.e. the a-h hne in Fig. ^6)). In this case, all the spectral densities are equal, Air^{Ljj) ~ A{lu), 
and can be calculated explicitly: 

1 r/2 
" ^ - + E.(r./2){i + iN- i)fAco)W ■ 

Here N denotes the degeneracy of the resonance and the level renormalization e{uj) = Eq + 
{N - l)J2ri'^r/2Tr)lTi[Ec/ma.-x{2TrkBT,eV/2,\u}\}]. For a single level without spin, N = 1, 
only charge fluctuations are possible and Eq. (|^) reduces to the Breit-Wigner formula with 
an unrenormalized Lorentzian form. For two spinless levels, N = 2, and the orbital index 
i behaves like a pseudospin; for two spin-degenerate levels we have = 4. If iV > 2, the 
(pseudo)spin dynamics contributes to Eq. (|^) and, at low temperature and for small excitation 
energies, it dominates over the charge fluctuations. As a consequence an additional sharp 
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Fig. 2 - o) Differential conductance Gi through dot 1 as a function of the bias voltage V = Vi = V2 
and for a fixed eia ~ — 4r. Temperature ksT = r/50. The two lower curves are the cotunnelling 
contributions to the conductance for Ae = (dotted) and Ae — T/2 (dashed). The higher curves 
are results in the resonant-tunnelling approximation the same sets of parameters (dot-dashed and 
solid, respectively) displaying zero-bias maxima as well as two satellites for Ae ^ 0. b) and c) show 
Aicr{uj) between two tunnel barriers; also the chemical potentials at Uri = ±eVi/2 are indicated. In 
6) Ae = and Vi = V2 = V, while in c) Ae / (still Vi = V2). The four solid Unes denote the 
energies of the states 1 |, 1 J,, 2 |, and 2 |. The states at eio- give rise to the broad Lorentzian shape 
in Aic,{u}). The positions of the sharp Kondo resonances in Aia{iL>) are depicted as the dashed lines 
(dotted lines denote the case Vi = 0), and the other relevant energies are shown on the sides. 

resonance, the Kondo peak, emerges in A{uj) at the positions of the chemical potentials in the 
leads. Figure |^&) shows an example of A{ui) for iV = 4 and introduces a schematic way for 
visualizing its relevant structure. The Kondo peak is reflected in the conductance as a sharp 
zero-bias maximum, dot-dashed curve in Fig. j^a), which is a distinctive sign of the Kondo 
effect. The peak in G is the higher the larger N is, in accordance with the A^-dependent 
Kondo temperature, see below. 

The most intriguing results are found for the general case with non-degenerate levels 
(£2(T ^ £i(T for detuned dots, and/or ^ Sii in a magnetic field) and/or different bias 
voltages, Vi ^Vi. In order to account for these extensions, we calculate Air,(uj) numerically. 
For brevity of the following discussion, we set the magnetic field to zero and focus on the 
conductance G\ through dot 1. Figure ||a) (solid line) shows an example of G\ for equal 
bias voltages, Vi = V, and a finite level separation Ae. The conductance reflects the shape 
of the spectral function ^1^(0;), illustrated in Fig. ^c), which exhibits resonances at ±eVi/2 
and — Ae ± eV2/2. The former arise due to the spin fluctuations in dot 1 and are similar to 
the usual Kondo peaks; the latter are due to tunnelling processes involving both levels - i.e. 
fluctuations of the orbital index i - and thus depend on the level separation Ae. As a rule, 
the conductance Gi shows a peak whenever two of the resonances in Aia-{uj) coincide: there 
is a zero-bias maximum due to the spin-degeneracy of the levels (although the levels are now 
detuned) and, in addition, two satellites at eV = ±Ae - these correspond to the additional 
peaks in Ai„{uj) and reflect the presence of the other level/dot. 

An experimental observation of the satellites in Gi is facilitated by two points. First, 
finite bias voltages induce decoherence to the electron transport and tend to suppress the 
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Fig. 3 - a) Linear conductance gi as a function of eV2 for eia- = — 4r and £20- = — 4.25r (dotted), 
-4r (solid), -3.94r (dashed) Q, -3.75r (long dashed), and -3.5r (dot-dashed), b) Differential 
conductance Gi at ksT = F/IO and the levels tuned to the mixed-valence regime - the renormalized 
upper level crosses the Fermi energies. The bare levels are at eia ~ — S.SF (solid), — 3.375r (dotted), 
and — 3.25F (dashed) and Ae = F; the renormalization is roughly 2.4F upwards. The energy scales F 
and ksT are also shown. 

Kondo effect for too high voltages However, in the present case the splitting Ae can 
be controlled starting from zero and only small bias voltages are required in observing the 
splitting. Second, due to the larger number of states (=4) the split peaks are more pronounced 
in the present system than, e.g., for a single dot with one Zeeman-split level (altough they 
have been observed 0] in this case as well). 

The fact that the voltages Vi can be chosen independently opens new interesting possibili- 
ties. Let us consider the spectral density in Fig. ||6) and the splitting of the Kondo resonance 
by the applied bias voltage. In the case of single quantum dots, the voltage-induced split in 
A{lu) is reflected as a mere monotonous decrease of the differential conductance G. In the 
present system, however, the linear conductance gi is found to show a peaked structure as a 
function of the voltages V2 and Vg2 (~ £20-) applied to the other dot, see Fig. ||a). The peaks 
in gi occur when the condition eV2 — ±2Ae is fulfilled, see Fig. ||c) and above. 

Before embarking on the scaling analysis, let us make a remark concerning the "Kondo 
regime" : it turns out that Gi may display a peaked structure also outside this regime. Let 
us consider a higher temperature than above, e.g. k^T = F/10, and tunc the system into 
the mixed-valence regime with the levels brought into a close vicinity of the Fermi energies 
in the leads. In this case, no Kondo effect is expected. A perturbative calculation produces 
the usual Coulomb blockade peaks, which in this case are four peaks that shift with eia\ the 
non-perturbative result for Gi shown in Fig. ^&) displays a triple-peak structure with one 
peak at zero-bias and two satellites at eV — ±Ae. The number and position of these peaks 
suggest that the peaks in Gi could be due to the Kondo effect after all, cf. results obtained 
above. On the other hand, if the peaks were of the Kondo origin, their width should reflect 
maxIfceT, fceTk} (or r~^, the inverse decoherence time of Ref. should this exceed the 
temperature), where Tk is the new low-energy scale, the Kondo temperature. However, in 
Fig. |6) the resonance widths are seen to be almost F, see the energy scales indicated in 
the figure; see also Rcf. Hence, although the peak widths are not a strict measure of 

the Kondo regime, we conclude that higher-order tunnelling processes may also give rise to 
non-trivial peaked structure other than the Kondo type. 

Scaling. In this section, we complement the above considerations in the framework of 
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the poor man's scaling technique [ p^ , a renormahzation group approach suited for the Kendo 
regime at equilibrium (zero-bias limit). This approach is used to derive explicit expressions 
for the logarithmic temperature dependence of the linear conductance g - experimentally the 
most convincing sign of the Kondo effect - and the Kondo temperature Tk. 

In order to obtain the proper scaling form of the conductance, we consider the dots deep in 
the Coulomb blockade regime, Sia <^ ^ri = 0, and start by Schrieffer- Wolff (SW) transforming 
the Hamiltonian Ti into the scattering form 



r.r' k,k' iaja' 



jicr.ja' t t }_ Tif'jf' t t 

'^rk,r'k'^rki<j'^r'k' ja'^ja''^ia- ~ J^'^rk.r'k' ^rkja' ^r' k' ja' '^irr'^iiT 



(4) 



Ticrjcr' _ \rpka (rpk' a' \* ( ^ \ I \ \_\ (e\ 

'^rk,r'k' ~ o ri r' j ) [ ^ , ^ ^ , p ^ I ■ W 

For N degenerate states, Ae = 0, Tiscat corresponds to the Coqblin-Schrieffer model [Q. The 
current through both dots, I = Ii + 12, calculated in the lowest (second) order in J is 

ia ja' 

where are the same probabilities that entered Eq. and 

ari-^r'j\^) — (■Jri.r'jPo) T^T ; Tj 7 {' ) 

exp[/3(cJ -Mr» +Mr'j)] - 1 

are the tunnelling rates from reservoir ri to r'j p^ . Here T^"^ ~ T^^i is assumed, implying 

•^r^r'j = Jrk'r'k'- should bc notcd that despite its different appearance, the present for- 
mulation yields just the cotunnelling current expressed above in terms of Eqs. (|l|) and (^. 
However, Eqs. (^ and (|^) are more suitable for the scaling analysis to follow. 

In the next step towards g{T), we perform poor man's scaling to second order in the 
coupling constants. In the degenerate case, Ae = 0, there is just one coupling constant J in 
the problem and it is renormahzed as J = J[l + NpoJlii^]~^ = [Npoluj^p^]^^ jl^]. The 
initial high-energy cutoff D corresponds to the band width in the leads, while the final cutoff 
is taken as D = max{kBT, eV} (or D = t~^, the decoherence rate of Ref. [^). The second 
equality for J follows from the definition of the Kondo temperature k^T^^ = D eyi'p[—l / N p^J]. 
This in turn yields an increasing Tk for increasing N and gives a qualitative reason, why the 
zero-bias peak in Gi is emphasized for larger values of iV, see above. 

Assuming that the zero-bias limit can be modelled in terms of equilibrium properties of 
the system, we insert the scaled coupling constants into Eq. (|^) and obtain 

for the leading cutoff dependence of the peak conductance. The conductance gi{T) through 
just one dot is sensitive to the voltage V2: If Vi — V2 — V, gi{T) = g{T)/2, while for V2 = 0, 
gi{T) is obtained from Eq. (^) by replacing the TV-dependent part by (N"^ + 4N - 2)/4iV^. 

Summary. In conclusion, we have studied electron transport through a capacitively- 
coupled double quantum dot with a large level spacing and intra and interdot charging ener- 
gies. The system is found to exhibit the Kondo effect due to fluctuations in both the orbital 
and spin degrees of freedom. The linear and nonlinear conductance is calculated for arbitrary 
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level positions and bias voltages and it is found to display a rich structure both in and out 
of the Kondo regime. For detuned levels, the differential conductance displays zero-bias as 
well as satellite peaks. The side peaks are argued to be experimentally observable despite the 
suppressive effect of finite bias voltages Q . The Kondo effect in each dot is shown to be very 
sensitive to the voltages applied to the other dot allowing delicate manipulation of the Kondo 
physics. An example of voltage-induced splitting of the Kondo peak has been discussed. 

* * * 

We would like to acknowledge J. Schmid, J. Weis, and U. Wilhelm for discussions. This 
work has been supported by the Finnish Cultural Foundation, EU TMR network "Dynamics 
of Nanostructures" , the Swiss National Foundation, and DFG through SFB 195. 

REFERENCES 

[1] H. Schoeller, in Mesoscopic Electron Transport, NATO ASI, Series E, Vol. 345, Kluwer (1997), 
and references therein. 

[2] T.K. Ng and RS. Lee, Rhys. Rev. Lett., 61 (1988) 1768; L.I. Glazman and M.E. Raikh, JETP 

Letters, 47 (1988) 452. 
[3] Y. Meir, N.S. Wingreen and R.A. Lee, Rhys. Rev. Lett., 70 (1993) 2601. 

[4] J. Konig, H. Schoeller and G. Schon, Rhys. Rev. Lett., 76 (1996) 1715; J. Konig et ai, Rhys. 

Rev. B, 54 (1996) 16820. 
[5] T. Rohjola et al, Europhys. Lett., 40 (1997) 189; J. Low Temp. Rhys., 118 (2000) 391. 
[6] W. Izumida, O. Sakai and Y. Shimizu, J. Rhys. Soc. Jpn., 67 (1998) 2444; A.L. Yeyati, F. Flores 
and A. Martin-Rodero, Rhys. Rev. Lett., 83 (1999) 600; M. Rustilnik, Y. Avishai and K. Kikoin, 
Rhys. Rev. Lett., 84 (2000) 1756. 
[7] D. Goldhaber-Gordon et al, Nature (London), 391 (1998) 156; S.M. Cronenwett, T.H. Oost- 

erkamp and L. Kouwenhoven, Science, 281 (1998) 540. 
[8] U. Wilhelm and J. Weis, Rhysica E, 6 (2000) 668; unpublished. 
[9] A. Kaminski, Yu.V. Nazarov and L.I. Glazman, Rhys. Rev. Lett., 83 (1999) 384. 
[10] J.M. Elzerman et al, J. Low Temp. Rhys., 118 (2000) 375. 
[11] H. Schoeller and J. Konig, Rhys. Rev. Lett., 84 (2000) 3686. 
[12] These probabilities can be obtained, e.g., by using detailed-balance relations. 
[13] The peak is maximized for a finite level splitting, e\a < E2a , due to the increase in the occupation 

probability of dot 1 and the subsequent increase in the weight under A\a{s^). 
[14] A.C. Hewson, The Kondo Problem to Heavy Fermions, Cambridge University Rress, 1993. 
[15] The functional form of these rates is similar to that for a single tunnelling barrier between two 

metal leads, emphasizing the idea of fluctuations extending through the dot. 
[16] In the non-degenerate case, Ae 7^ 0, the scaling equations would need to be solved numerically. 
This is also the case close to the point h of Fig. nl6); there the coupling constants are found to 
scale, but in a different way. 



